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Abstract 

03 _ 

We give sharp bounds for the isotropic unimodal probability convolution semigroups 
when their Levy-Khintchine exponent has weak local scaling at infinity of order strictly 
between and 2. 



^-i 1 Introduction 

Estimating Markovian semigroups is important from the point of view of theory and applications 
because they describe evolution phenomena and underwrite various forms of calculus. Diffusion 
semigroups traditionally receive most attention [2A\ but considerable progress has also been 
made in studies of transition densities of rather general jump-type Markov processes. Such 
! studies are usually based on assumptions concerning the profile of the jump or Levy kernel 
(measure) at the diagonal (origin) and at infinity [TTJ [H]. As a rule the assumptions can be 
viewed as approximate or weak scaling conditions for the Levy density, to which some structure 
conditions may be added, see [Tlj (1.9)-(1.14) and Theorem 1.2]. Typical results consist of sharp 



in . 

O . two-sided estimates of the heat kernel for small and/or large times. 

Transition semigroups of Levy processes allow for a deeper insight and direct approach from 
several directions thanks to their convolutional structure and the available Fourier techniques. 
For instance, the upper bounds for transition densities of isotropic Levy densities with relatively 
fast decay at infinity are obtained in [201 [21] by using Fourier inversion, complex integration, 
saddle point approximation or the Davies' method. In this work we study one-dimensional 
distributions pt{dx) of rather general isotropic unimodal Levy processes X = (X t , t ^ 0) in 
M. d [26] . We focus on pure-jump isotropic unimodal Levy processes. Thus, X is a cadlag 
stochastic process with distribution P, such that X(0) = almost surely, the increments of X 
are independent with rotation invariant and radially nonincreasing density function pt(x) on 
M. d \ {0}, and the following Levy-Khintchine formula holds for £ e M. d , 

Ee i{$,xt) = f e i(M pt (dx) = e~ tm , where ^(0=/ (1 - cos (£,x))v{dx). (1) 
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Here and below v is an isotropic unimodal Levy measure and E is the integration with respect 
to P. Further notions and definitions are given in Sections [2] and [3] below. Put differently, we 
study the vaguely continuous spherically isotropic unimodal convolution semigroups (p t , t ^ 0) 
of probability measures on M. d with purely nonlocal generators. (In this work we never use 
probabilistic techniques beyond the level of one-dimensional distributions of X.) 

Our main result provide estimates for the tails of p t and its density function pt(x), expressed 
in terms of the Levy-Khintchine exponent ip. We also use ip to estimate the density function 
v(x) of the Levy measure v. Since ip is radially almost increasing, it is comparable with its 
radial nondecreasing majorant ?/>*, and as a rule we employ ip* in statements and proofs. The 
extensive use of ip (ip*) rather than v is a characteristic feature of our development and may be 
considered natural from the point of view of pseudo-differential operators and spectral theory 
[16J. As usual, the asymptotics of ip at infinity translates into the asymptotics of v and p t at 
the origin. Our estimates may be summarized as follows, 

p t (x)^[r(l/t)] d A tr ^P «ft(0)AH g )], (2) 

see Theorem [2H Corollary [23] and ff24|) for detailed statements. Here means that both sides 
are comparable i.e. their ratio is bounded between two positive constants, ip is assumed to 
satisfy the so-called weak upper and lower scalings of order strictly between and 2, and ip~ is 
the generalized inverse of ip*. The bound ([2]) holds locally in time and space, or even globally, 
if the scalings are global. We note that the corresponding estimates of u, to wit, 



u(x) ~ 



ip*(\x 


n 




x\ 


d 



(3) 



are simply obtained as a consequence of (j2J) for small time, and not as an element of its proof 
(see Corollary |2"3"]) . It is common for v to share the asymptotics with p t because v = lim t _>oPt) 
a vague limit in M. d \ {0}. It is also a manifestation of the general rule mentioned above that 
?/>*(|x| _1 ) in our estimates reflect the properties of pt(x) and v(x). The denominator, \x\ d , 
in the estimates comes from the homogeneity of the volume measure in M. d (see the proof of 
Corollary [23]) and ip~(l/t) approximates pt{0) as follows from a change of variables in Fourier 
inversion (cf. [26] Lemma [16] and Lemma ITT)) . All these reasons make the above bounds most 
natural. Therefore below we shall address (two-sided and one-sided) estimates similar to 02]) 
and (EJ) as common bounds. We note that the common upper bounds, 

ftfrKg*™, „(,)<cgflM, X e^\{0},t>0, 

hold with a constant depending only on the dimension for all isotropic unimodal Levy semi- 
groups. This is proved in Corollary [TJ The lower common bounds hold if and only if tp has (the 
so-called weak) lower and upper scalings of order strictly between and 2. Indeed, in Theo- 
rem [22] we show that for unimodal Levy processes, the scaling of ip (at infinity) is equivalent to 
the common bounds for the transition density and the Levy measure (at the origin). In fact, 
already the lower bound v(x) ^ c^*(|x| _1 )/|a;| d implies such scalings of ■0. 

We thus cover all the cases of isotropic Levy processes with scaling studied in literature, 
and upper bounds are provided for all isotropic unimodal Levy processes, which (only) leaves 
open the problem of estimating the semigroup of (scaling) unimodal Levy process with Levy- 
Khintchine exponent slowly varying at infinity or regularly 2- varying at infinity |6][22]. Here typ- 
ical examples are the geometric stable processes and the variance gamma process-the Brownian 
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motion subordinated by an independent gamma subordinator. As we see from [61 Section 5.3.4], 
such processes require specialized approach and their transition density and Levy-Khintchine 
exponent do not easily explain each other. 

Bochner's procedure of subordination is strongly rooted in semigroup and operator theory, 
harmonic analysis and probability [32j [27J [26] . In the present setting it yields a wide array 
of asymptotics of ip, pt and z/, which explains intense current developments. In particular, 
common bounds were recently obtained for a class of subordinate Brownian motions, mainly 
for the complete subordinate Brownian motions defined by a delicate structure condition [26J . 
Highly sophisticated current techniques and results in this direction are presented in |T8] , see 
also [6l [26]. Our approach is, however, more general and synthetic; we demonstrate that the 
sources of the asymptotics of (unimodal) Pt{x) are merely its radial monotonicity and the 
scalings of ijj, rather than further structure properties of ip. 

We illustrate our results with several classes of relevant examples. These include situations 
where former methods cannot be easily applied and the present method works well. In this 
connection we note that ip is an integral quantity and may exhibit less variability than z/, in 
particular the scaling properties of ip are more easily manageable. Many of our examples are 
in fact the subordinate Brownian motions, and then = 0(l£| 2 )> where <fi is the Laplace 
exponent of the subordinator, i.e. a Bernstein function. There is by now an impressive pool of 
Bernstein functions studied in the literature, with distinct asymptotics at infinity. For instance 
the monograph [2B] gives well over one hundred cases and classes of Bernstein functions in 
its closing list of examples. Many of these functions have the scaling properties used in our 
paper. This immediately yields sharp estimates of the Levy measure and transition density 
of subordinate Brownian motions corresponding to such subordinators. In comparison, former 
methods require to first find estimates of the Levy measure of the subordinator (this is where 
the completeness of the subordinator plays a role), then to estimate the Levy density of the 
corresponding subordinate Brownian motion and then, finally, to estimate its semigroup [TU| 

USE]. 

Our results and arguments are purely real-analytic. We circumvent those additional steps 
and directly estimate the heat kernel by arguments not unrelated to integration by parts. In 
particular, for subordinate Brownian motions with scaling, we relax the usual completeness 
assumption on the subordinator. Noteworthy, while on one hand the completeness need not be 
used when scaling is present, on the other hand the complete subordinate Brownian motions 
exhibit all the types of asymptotics of ip, pt and v of general unimodal Levy processes with 
scaling. This is in particular manifested in Corollary [271 

We remark that analogues of the on-diagonal term [ip~(l/t)] d w Pt(0) in ([2]) are often 
obtained for more general Markov processes via Nash and Sobolev inequalities [3j [8j [2j [291 EI] • 
For our unimodal Levy processes we instead use Fourier inversion and (weak) lower scaling, 
see Proposition [T9l Also our approach to the off-diagonal term t?/;*(|x| _1 )/|x| d is very different 
and much simpler than the arguments leading to the upper bounds in the otherwise more 
general Davies' method [8J, [121 Section 3]. Our common upper bounds are straightforward 
consequences of a specific quadratic parametrization of the tail function, which is crucial in 
applying the techniques of Laplace transform. The common lower bounds are harder, and they 
are intrinsically related to upper and lower scalings via certain differential inequalities in the 
proof of Theorem [26j 

The (local or global) comparability of the common lower and upper bounds is a remarkable 
feature of the class of semigroups captured by Theorem [221 We expect further applications of 
the estimates. For instance, under (weak) global scalings we obtain important metric-type [TTj 
global comparisons P2t(x) ~ Pt(x) andp 4 (2x) ~ Pt(x), given in Corollary 1241 b elow . These should 
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matter in perturbation theory of Levy generators and in nonlinear partial integro-differential 
equations. Since uniform estimates are important in some applications, in Corollary [2H and 
elsewhere in the paper the comparability constants are shown to depend in a rather explicit 
way on specific properties of the semigroups, chiefly on scaling. For instance for the isotropic 
a-stable Levy semigroup in M. d with < a < 2 |6J, we have ^(O = lC| a > an d we arrive at 

with explicit constants given by ( I17j) . ( 1261) and (1291) below. Noteworthy, our (weak) scaling 
conditions imply majorization and minorization of ip at infinity by power functions with ex- 
ponents strictly between and 2, but do not require its comparability with a power function, 
see examples in Section 14.11 Furthermore, the exponents a and a in the assumed scalings only 
affect the comparison constants in common bounds, but not the rate of asymptotic, which is 
solely determined by ip. 

For convolution semigroups of probability measures more general than unimodal, the struc- 
ture of the support and the regularity of the Levy measure plays a crucial role. In particular 
the directions which are not charged by the Levy measure see in general lighter asymptotics of 
Pt(dx) [l_3j [31]. In consequence, the estimates of severely anisotropic convolution semigroups 
require completely different assumptions, description and methods. Our experience indicates 
that v surpasses if) in such cases. Estimates and references to anisotropic v with prescribed 
radial decay and rough spherical marginals may be found in [28] (see also [31] for more details 
in the case of homogeneous anisotropic v). 

The structure of the paper is as follows. In Section [2] we discuss first consequences of 
isotropy and radial monotonicity. In particular we compare ip with Pruitt-type function hi(r) = 
J Rd (1 A \xi\ 2 /r 2 ) v{dx) and we estimate from above the tail function f t (p) = F(\X t \ ^ yfp) by 
using Laplace transform and ip. This and radial monotonicity quickly lead to upper bounds for 
Pt{x) and v{x). In Section [3] we discuss almost monotone and general weakly scaling functions. 
In SectionH]we specialize to scalings with lower and upper exponents a, a G (0, 2), and we give 
examples of ip with such scaling. To obtain lower bounds for pt(x) and v{x), in Lemma [TBI we 
recall an observation due to M. Zahle, which is then used in Lemma [T41 to reverse the comparison 
between ip and the tail function f t . The generalized inverse ip~ plays a role through a change of 
variables in Fourier inversion formula for pt(0) in Lemma [161 and through equivalence relation 
defining "small times" (stated as Lemma [17]) . In Theorem I2T1 we combine all the threads to 
estimate p t , as summarized by (T5]). In Corollary [23] we obtain ([3]) as a simple consequence of 
Theorem [2TJ To close the circle of ideas, in Theorem [22] we show the equivalence of the weak 
scalings with the common form of bounds for p t and v. In Proposition [2S] we state a connection 
between v and ip for a class of approximately isotropic Levy densities. 



2 Unimodality 

We shall often use the gamma and incomplete gamma functions: 

/*oo pt poo 

T(S)= / e'V- 1 ^, t(M)= / e _ V -1 d«, T(5,t) = / e^u^du, 5, t > 0. 
Jo Jo Jt 

Let M. d be the Euclidean space of (arbitrary) dimension d G N. For x G M. d and r > we let 
B(x,r) = {y G R d : \x - y\ < r}, and B r = 5(0, r). We denote by oo d = 2n n ' 2 /T(n/2) the 
surface measure of the unit sphere in R d . All sets, functions and measures considered below are 
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(assumed) Borel. A (Borel) measure on M. d is called isotropic unimodal, in short: unimodal, if 
on M d \ {0} it is absolutely continuous with respect to the Lebesgue measure, and has a finite 
radial nonincreasing density function. Such measures may have an atom at the origin: they 
are of the form a5o(dx) + f(x)dx, where a ^ 0, 5q is the Dirac measure, 

POO 

f(x)= l Br (x)ii(dr) = ii((\x\,oq)) (a.e.), 



and fJ,((£, oo)) < oo for all e > 0. A Levy process X = (X t , t ^ 0), is called isotropic unimodal 
(in short, unimodal) if all of its one-dimensional distributions (transition densities) pt{dx) are 
such. Recall that Levy measure is any measure concentrated on R d \ {0} such that 

(\x\ 2 A l) vidx) < oo. (5) 



Unimodal pure-jump Levy processes are characterized in [30J by unimodal Levy measures 
u(dx) = v[x)dx = u{\x\)dx. 

Unless explicitly stated otherwise, in what follows we assume that A is a pure-jump uni- 
modal Levy process in W 1 with (unimodal) nonzero Levy measure (density) v. 

Each measure pt is the weak limit of 



ft = 



oo 1 
n=0 



where e — > + and z/* n are convolution powers of measures v e (A) = u(A \ B e ), see, e.g., [6j 
Section 1.1.2]. Each p t has a radial nonincreasing density function pt{x) on M. d \ {0} and atom 
exp[— tu(M. d )] at the origin if u(M. d ) < oo (no atom if tv is infinite). 
For r > we define after [23] . 

h(r) = / A X ) ' L ( r ) = v ( B r) ■ ( fi ) 

Clearly, ^ L(r) < h(r) < oo, L is nonincreasing and h is decreasing. The strict monotonicity 
and positivity of h follows since v ^ is nonincreasing, hence positive near the origin. 

The first coordinate process X\ of X t is unimodal in R. The corresponding quantities L\ (r) 
and hi{r) are given by the (pushforward) Levy measure V\ = uon^ 1 , where n 1 is the projection: 
M. d 3x = (xi, . . . , Xd) xx, see [251 Proposition 11.10]. With a typical abuse of notation we let 
vi{y) denote the (symmetric and nonincreasing on (0, oo)) density function of V\\ 



My) 

Thus, 



/ KvV + \A 2 )dz, yGM\{0}. 



//.,.(/•) = / ( 4 A A Mdy) = l (^4- A 1 ) "(^)- r > 0. 



q^2i J J \ ^-i2 



Therefore, 

hi(r) ^ h(r) ^ h(r)d, r > 0. (7) 
In fact, (J7|) is valid more generally: for all rotation invariant Levy measures. 
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Remark 1. h(r) and h\(r) are decreasing, while r 2 h(r) and r 2 h\(r) are nondecreasing. 

Since ^ is a radial function, we shall often write ip(u) = where (6R d and u = |£| ^ 0. 
We obtain the same function for X}. Clearly, tp(0) = and, as before for h, ip(u) > for u > 0. 
We now show how to use hi to estimate the Levy-Khintchine exponent ip of X. 

Lemma 1. 

l hl ( I ) ^ <: 2h x (-), u> 0. (8) 



TT 2 \ U I \U 



Proof. For t ^ we define n(t) = J *(l — cos r)dr, and we claim that 

(r 2 A l)dr sC < 2 / (r 2 A l)dr, t ^ 0. (9) 



7T 2 



Indeed, 1 — cosr = 2sin 2 (r/2) ^ 2(r 2 A 1), which gives the upper bound. If ^ r ^ t ^ n, 
then 1 — cosr = 2sin 2 (r/2) ^ 2r 2 /n 2 hence, n(t) ^ 27r~ 2 J '(r 2 A l)<ir, and if t > 7T, then 



/ (1 - cos r)dr = (i - tt/2) + (1 - sin t) ^ f dr, 

J-k/2 Jtt/2 



which yields (Q. (The constants in ([9]) may be improved.) We define an auxiliary measure \i 
on (0, oo) by letting 

oo)) = Mv) for a - e - y > °- 

Let w > 0. By a change of variables and Fubini-Tonelli, 



1, (\ 



oo 

( - ) = / [(A 2 ) A \\v x {y)dy = I \{u 2 y 2 ) A 1] / fx(dt)dy 



OO /*f /* OO "1 pt/'UL 

J [(u 2 y 2 ) A l]dy fi(dt) = / -/ [r 2 A l]<ir fi(dt). 
o Jo Jo u Jo 

Similarly, 

I poo poo p poo 

= / [1 - cos(«y)] v\{y)dy = / [1 - cos(«j/)] / fjt(dt)dy = / -K(tu)/i(dt). 

* Jo Jo J(v,oo) Jo u 

By these identities and (jHJ) we obtain (jSJ). □ 

We define the maximal characteristic function if)*(u) : = sup ss ; u i>(s), where u ^ 0. The 
following result is a version of [El Proposition 1]. 

Proposition 2. 

^(it) < i/)*(u) ^ 7r 2 ^(n) /or u ^ 0. (10) 
Proof. Since is nonincreasing, by Lemma [1] for u ^ we have 

«S ip*(u) < 2 sup /ii f - J = 2/ij ( - J < tt 2 ^(u). 
o<s<u V s / 

□ 
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We write f(x) ~ g(x) and say / and g are comparable if /, g ^ and there is a positive 
number C, called comparability constant, such that C _1 /(x) ^ g(x) ^ C f(x) for all x. We 
write C = C(a, . . . , z) to indicate that C may be so chosen to depend only on a, . . . ,z. We 
say the comparison is absolute if the constant is absolute. Noteworthy, while tp is comparable 
to a nondecreasing function, it need not be nondecreasing itself. For instance, if ip{u) = 
[1 - (sinu)/u], then ^'(2tt) = -| < 0. 

The following conclusion may be interpreted as relation of "scale" and "frequency". 

Corollary 3. We have h(r) &hi(r) « ip(l/r) m ip*(l/r) for r > 0. 

Proof. The constant in the leftmost comparison depends only on the dimension, see flTJ). The 
other comparisons are absolute, by Lemma [T] and Proposition HJ □ 

By Corollary [3] and definitions of L%, L and h, we obtain the following inequality, 

Li(r) < L(r) < h(r) ^ Cty*(l/r), r > 0. (11) 

Our main goal is to describe asymptotics of u(x) and pt{x) in terms of ip*. We start with 
analysis of the Laplace transform of the integral tails of pt- For reasons which shall become 
clear in the proof of the next result, we choose the following parametrization of the tails, 

f t (p)=n\X t \ > v^) = P(|X t | 2 >p), p^Q, t>0. (12) 

We consider the Laplace transform of f t , 

POO 

Cf t (X) = / e- x ?f t (p)dp, A ^ 0. 
Jo 

Lemma 4. There is a constant C\ = C\{d) such that 

(l - e-^)) ^ Cf t (X) ^ (l - e*TW) , A > 0. 

Proof. By Fubini's theorem, J Rd h(x)k(x)dx = L d h(x)k{x)dx for integrable functions h, k. By 
this, (fT2l) and change of variables we obtain, 

A£/,(A)=E(l-e- A l x 'l 2 ) 

= 1- / e- x W°pt(x)dx = (4vr)- d / 2 / fl - e "^M) e~^dx. 

JR d JR d ^ ' 

By [T5l Theorem 2.7], 

ip(su) ^ ip*(su) ^ 2(s 2 + 1)V>*H s, u ^ 0. (13) 
(The estimate may usually be improved for specific if).) We also note that 

l-e- bt <&(l-e -t ), t^0, b^l, (14) 

and we obtain, 

XCf t (X) ^ (4tt)^/ 2 [ (l- e -*(M a +Wfr/XA e -l-l 2 /4 rfx 

^ 2(2d+l) fl -e-^'^V (15) 
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On the other hand, if |x| ^ 1, then ip (xVtj ^ ip* (\x\Vx) /tt 2 ^ ip* (ytj /tt 2 by flE]). Thus, 
XCft{\) > (4vr)- d / 2 [ e-\*\ 2 /*dx(l- e^W^) > (l- e"*^ , 

y B e v / 7r 2 r(d/2) v / 

where we use (JT3J) and the upper incomplete gamma function ■). □ 

The upper bounds for tails shall follow from this auxiliary lemma. 
Lemma 5. If f is nonnegative and nonincreasing, then for n, m = 0, 1, 2, . . . and r > 0, 

< —-—!——— r~ n-m '~ 1 1 (£[s"Y]) (ri) (r _1 ) | . 
7(n + m + 1, 1) 

Proof. If U > and r = u^ 1 , then 

u n+m+1 |(£[s"7]) (n) (u)| = u n+m+1 |(-l) n £[s n+ "7](u)| 

^ u n+m+l / e- s?t s n+m /(s)c/s ^ /(w- 1 ) / M(Ms) n+m e- s?t rfs 
Jo Jo 

= /(m- 1 ) / M n+m e-"c/M = /(M- 1 ) 7 (n + m + 1,1), 



where we use the lower incomplete gamma function 7(-, •)• ^ 

The following estimate results from (fl5|) and Lemma |5] with n = m = 0. 

Corollary 6. For r > we have F(\X t \ ^ r) ^ J^L(2d+ 1) (l - e -*^ 1 / r >). 

Here is a general upper bound for density of unimodal Levy process. (As we shall see in 
Theorem [2T1 and |26| a reverse inequality often holds, too.) 

Corollary 7. There is C = C(d) such that p t (x) ^ Ct^*{l/\x\) /\x\ d for x G R d \ {0}. 
Proof. By radial monotonicity of y y-y p t (y), ffT3j) and Corollary [61 

! - r,- d n / lvl . . rf M*(l/|x|) 



P(|x 


1/2 <l 




< 


N) 




#|x| \ 


/2 





1 - 2" d )w d V 2 y 1 1 " \x 



\d 



□ 



Since p t (x)/t — y v(x) vaguely on R \ {0}, we also obtain 



v{x)^C r ^f\ i£K d \{0}. (16) 
Tracking constants, e.g., for the isotropic a-stable Levy process addressed in (j3J), we get 

d^a-p (d+a) t 
Pt(x) < —, ., V .? '-. r ^J~, t > 0, X E R d . (17) 

F v ' 2(1 - 2- d )7r d / 2 7 (1, 1) \x\ d+al ' v ' 

(In fact, for this constant we override (fl3|) by ip(su) = s a ip{u) in the proof of Lemma ffl) 



<S 



3 Weak scaling and monotonicity 

For reader's convenience we shall give a short survey of almost monotone and weakly scaling 
functions. The results may be considered folklore [18] , but the actual variants which we need 
may be difficult to find. 

Let <p : / — > [0, <x>], for a connected set I C [—00,00]. First, we call <p almost increasing if 
there is (oscillation factor) c G (0, 1] such that c(f>(x) ^ <f)(y) for x,y G I, x ^ y. Let 

4>*{y) = snp{4>(x) : x G /, x ^ y}, y G /. 

We easily check that <fi* is nondecreasing, <fi ^ <p* and the following result holds. 

Lemma 8. <ft ^ s almost increasing with oscillation factor c if and only if c<p* ^ (f>. 

E.g. the Levy-Khintchine exponent if) in Proposition [2] is almost increasing with factor l/V 2 . 
On the other hand, if there is C G [1, 00) such that C(f>(x) ^ (f>(y) for x,y G / , x ^ y, then we 
call (ft almost decreasing (with oscillation factor C). Let 

<p*(x) = sup{<%) : y G /, y > x}, x el. 

We easily check that 0* is nonincreasing, 0^0* and the following result holds. 

Lemma 9. <ft ^ s almost decreasing with oscillation factor C if and only if <p* ^ C(f). 

We note that <fi is almost increasing on / with factor c if and only if is almost decreasing 
on / with factor 1/c. Here is another simple observation which we give without proof. 

Lemma 10. Assume that sets ii, J2, / — I\ U I2 C M. are connected. If <fi is almost increasing 
(decreasing) on I\ with factor d (C), almost increasing (decreasing) on I2 with factor c" (C"), 
then <p is almost increasing (decreasing) on I with factor c = c'c" (C = C'C" '). 

We now specialize to / = (9, 00), where 9 G [0, 00). We say that (f> satisfies the weak lower 
scaling condition (at infinity) if there are numbers aG IR. and cG (0, 1], such that 

<j}{\9) ^ c\^(f)(9) for A ^ 1, 9e I. (18) 

In short we say that (f> satisfies WLSC(a, 9, c) and write (f> GWLSC(a, 9, c). If G WLSC(a, 0, c), 
then we say that <fi satisfies the global weak lower scaling condition. 

Similarly, we consider / = (9, 00), where 9 G [0, 00). The weak upper scaling condition 
holds if there are numbers «6 K. and CG [1, 00) such that 

(j){\9) ^ 0X^(^(9) for A ^ 1, 9e I. (19) 

In short, <fi G WUSC(a, 9, C). For global weak upper scaling we require 9 = in ( TT9l) . 
The following clarification is an analogue of (5j Theorem 2.2.2]. 

Lemma 11. We have <p G WLSC(a,9_,c) if and only if (f)(9) = n(9)9- and k is almost increasing 
on (9,oo) with oscillation factor c. Similarly, (f> eWUSC(a,9 ,C ) if and only if (f)(9) = n(9)9 a 
and k is almost decreasing on (9, 00) with oscillation factor C . 
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Proof. Let G [0, oo), I = (0, oo) and GWLSC(a,0,c). Let «(0) = 0(0)0~ a on J. If £ < 77 ^ 
and A = 0/77, then 

«(0) = 0(Ar/)(Ar/)^ ^ cA^0(r/)(Ar/)^ = ck(t]). 

Thus, k is almost increasing. On the other hand, if (f)(6) = k{9)9— and k is almost increasing 
with factor c, then for A ^ 1 and G I we have 

0(A0) = k(A0)(A0)^ ^ ck{9)(\9)^ = cA^0(0). 

Thus, GWLSC(a,0,c). The proof of the second part of the statement is left to the reader. □ 

In particular by Lemma [TIJ WLSC(O,0, c) characterizes almost increasing functions on 
(9, 00), and WUSC(0, 9, C) characterizes almost decreasing functions on (9, 00). 
For example, h,h x G WUSC(0,0,l)nWLSC(-2,0,l), see Remark[U 

Remark 2. If_0_GWLSC(a,0,c) and c0 ^ ip ^ C<f>, then^ GWLSC(a ,9, cc/C). Similarly, if 
GWUSC(a,0,C) and c0 < ip < C0, then tp GWLSC(a,0,C C/c). 

As or 9 decrease, the scaling conditions tighten. Here is a loosening observation. 

Lemma 12. Let : [9 1,00) — > (0, 00) be nondecreasing and nonzero. If < 9\ < 9 and G 
WLSC(a, 0, c), then G WLSC(a, 9 X , c x ) with c x = c (9 1 /9) 1 - 1 0(0 1 )/0(0) . If rather < 9 X < 9 
and G WUSC(a, 9,C), then G WUSC(a, 9 X , C x ) with C X = C (9/9 x ) 0(0)/0(0x) . 

Proof. In view of Lemma [TU1 and Lemma [TT| it is enough to study k(9) = 0(0)0~- on [9 X , 0]. If 
X ^ 9 ^ 0, then 0(0! ) A 0^-) k(0) < 0(0) (0^- V This gives the first implication 

and the second obtains similarly. □ 

Remark 3. Let ^ be continuous an increase to infinity. If G WLSC(a, 0, c) [WUSC(a, 0, C)], 

then0- 1 G WUSC(l/a, 0(0), c" 1 ^) [WLSC(l/a,0(0),C _1/ °), resp.]. Indeed, since is increas- 
ing, by the scaling for A ^ 1 and > we have 

A0~ 1 (0(0)) = A0 = 0- x (0(A0)) > 0~ 1 (cA^0(0)). 

Thus for (arbitrary) y = 0(0) > 0(0), if s — cA— ^ c, in particular if s ^ 1, then 

tr 1 '*f*<r 1 (v)><i>- l {sv)- 

Similarly, if G WUSC(a,0,C), then 

C^V/^- 1 ^) < 0- 1 ( S2/ ), y > 0(0), s > C. 
For 1 ^ s < C, by monotonicity of -1 , 

C-^s 1 ^- 1 ^) < 0- 1 (y) < 0- 1 (si/), y > 0(0). 

This proves our claim. 

Remark 4. We also note that GWLSC(a,0,c) if and only if 1/0 GWUSC(-a,0, 1/c). Sim- 
ilarly, 0(t) GWLSC(a, 0,c) if and only if 0(l/t) GWUSC(-a, 0, 1/c). For instance, our Levy- 
Khintchine exponent ip has global upper scaling: by Remark [TJ Lemma [11] and Lemma [TJ 
V>, V* gWUSC(2,0,tt 2 ), and ij> GWLSC(0,0,l/vr 2 ). In fact, by flU we have 

i)*{u) ^ ip*(\u) ^4AV(«), \>1,u^0, 
and so -0* GWUSC(2,0,4). Of course, -0* GWLSC(0,0,1), meaning that -0* is nondecreasing. 
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4 Scaling of the Levy-Khintchine exponent 



We shall study consequences of scaling of the Levy-Khintchine exponent if> of (isotropic) uni- 
modal pure-jump Levy process X with nonzero Levy measure v. For economy of notation, in 
the sequel we only consider (assume) scaling exponents a, a satisfying: 

< a < 2 and < a < 2. (20) 

4.1 Examples 

The Levy-Khintchine (characteristic) exponents of unimodal convolution semigroups which we 
present in this section all have lower or upper scaling suggested by (12"U|) . This can be verified 
in each case by using Lemma [TTJ While discussing the exponents, we shall also make connec- 
tion to subordinators, special Bernstein functions and complete Bernstein functions, because 
they are intensely used in recent study of subordinate Brownian motions, a wide and diverse 
family of unimodal Levy processes cf. [18J. The reader may find definitions and comprehensive 
information on these functions in [26]. When discussing subordinators we usually let </?(A) de- 
note their Laplace exponent, and then if)(x) = ip(\x\ 2 ) is the Levy-Khintchine exponent of the 
corresponding subordinate Brownian motion. We focus on scaling properties of if). 

1. Let if){u) = ^YTi=\ c i uai i where u ^ 0, n G N, Q > 0, G (0,2). This if) is the Levy- 
Khintchine exponent of X = Y^h=i ^"^(O' wnere Xu\ are independent isotropic te- 
stable Levy processes (in R d ). Let a = minaj, a = maxa^. Then for every e > 0, 
if) GWLSC(a, e, c), where c = Cj/(Y^ =1 CiE ai ~ a ) and j is such that a = atj. Even simpler, 
if) GWUSC(a,0, 1) and if) GWLSC(a,0, 1). The latter is also true if if)(u) = J u a fx(da), 
where \i is a finite measure on [a, a] (and < a < a < 2). More generally, let A be a 
(Borel) set and let ip(s,u) = ip s (u) satisfy WLSC(a,fl,c) [WUSC(a, 6, C)] for s G A, and 

MZ)=MS\)= f (l-coa(Z,x))v 8 (dx), £eR d , 

for some unimodal Levy measures v s (to wit, u.(-) is assumed jointly Borel). Let /i be a 
measure on A such that f. ip s (l)fj,(ds) < oo. Then ip(x) = f A ip s (x)fi(ds), which is finite 
by Remark |2j is the Levy-Khintchine exponent of a unimodal Levy process and we have 
if) G WLSC(a,^,c) [WUSC(a,0,C7), respectively]. 

2. Let ip(\) = / °°(1 — e~ Xu )n(dr) be a Bernstein function [26], i.e. the Laplace exponent of 
a subordinator r\ [261 HI 1251 EE] , and let Y be an independent (isotropic) unimodal Levy 
process with characteristic exponent x- Then the process X t = Y Vt is unimodal and 
has the the characteristic exponent if)(x) = (p(x( x )) [26J. If x WUSC(<5i, 9, C\) and 
if GWUSC(a 2 , x{@)i C2); then if> GWUSC(aia 2 , 9, C° 2 C 2 ). From concavity of Bernstein 
functions it also follows that if % G WLSC^, 0, c 1 ), = ini o^gx(9) and GWLSC(a 2 , 9*, c 2 ), 
then ^ gWLSC^^) £lC 2 ). We always have ^ x(0/tt 2 , see Proposition El and of- 
ten = Of particular interest here is x(C) — l£| 2 ; i- e - — B%t, where B is the 
Brownian motion in M. d . The process X is then called a subordinate Brownian motion. 
Furthermore, it is called special subordinate Brownian motion if the subordinator is spe- 
cial (i.e. given by a special Bernstein function), and it is called complete subordinate 
Brownian motion if the subordinator is complete [26]. The (unimodal) Levy measure 
density of Y is given by the formula 

v(x) = / {A7rty d / 2 e~-^fi{dt), (21) 
Jo 
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and its Levy-Khintchine exponent ^(|£| 2 ) is in WUSC(2a 1; 2 , c\) or WLSC(2a 1 , 6 2 , cf), 
respectively. 

3. Let ^(0 = |£| a log^l + |£| 7 ), where 7, a, a + 2(3 G_(0, 2). If e > 0, then if) G WUSC(a + 
e, 1, C) n WLSC(a - e, 1, c) for some < c < 1 ^ C < 00. Furthermore, if) e WUSC(a + 
7/3 + ,0, 1) fl WLSC(a— 7/3_, 0, 1). We note that if) is the Levy-Khintchine exponent of 
a subordinate Brownian motion, see Theorem 12.14, Proposition 7.10, Proposition 7.1, 
Corollary 7.9, Section 13 and examples 1 and 26 from Section 15.2 in [26]. Many more 
examples related to subordinate Brownian motions readily follow from [261 Section 15]. 

4. Let /i(r) = (1 — sinr)r _1_a / 2 , where a G (0,2) and — J °°(^ ~ e~ r ^ 2 )fi(r)dr; in fact, 
let f(t) be a nonnegative bounded function on (0, 00) and define /x(r) = f{r : )r~ 1 ~ a l 2 . Then 
ip(£) = Jq°°(1 — e~ r '^' )/j,(r)dr is the characteristic exponent of a subordinate Brownian 
motion and ip(0 ^ C|£| Q - We also note that 

poo r\£\~ 2 poo 

iP(£)= / (l-e^'Xrldr^e" 1 ^ r^(r)dr + (I - e' 1 ) fi(r)dr. 

Jo Jo J\£\- 2 

If liminf| €Koo |^| 2_Q J^' -2 f{r)r~ a / 2 dr > 0, then ^(0 > c\£\ a , hence ~ for 

|£| 1. For example f(r) = (1 — sinr) and f(r) = sin 2 (r + r -1 ) produce such if). For 
later use we note that these if) may be considered Levy-Khintchine exponents of subordi- 
nate Brownian motion but the Laplace exponents of the respective subordinators are not 
complete Bernstein functions, as their Levy measure densities //(r) are not monotone, cf. 
[261 Definition 6.1]. 

5. Let X be pure-jump unimodal with infinite Levy measure and Levy-Khintchine exponent 
if). Let a scaling condition with exponent a or a hold for if). For fixed r > 0, we let X r 
be the (truncated) unimodal Levy process obtained by multiplying the Levy measure of 
X by the indicator function of the ball B r , and let if) r be its Levy-Khintchine exponent. 
Since ^ if) — if) r is bounded, if> r is comparable with if) at infinity, and so if> r has (local) 
scaling with the same exponent as if). For later discussion we observe that if) r is not an 
exponent of a subordinate Brownian motion because the support of its Levy measure is 
bounded [26], Proposition 10.16]. 

6. We may combine the two preceding examples to obtain subordinate Brownian motions 
which have local scaling but are not special, since the Levy measure of their subordinators 
is being truncated. In fact, let a G (0,2) and n(dr) = ^2^ =2 ^i/k(dr) (k a l 2 — (k — l) a ^ 2 ), 
<p(X) = Jo (1 ~~ e~ rX )fi{dr). This ip is not a complete or even special Bernstein function 
[2"6j Proposition 10.16]. As usual, if>(£) = y?(|£| 2 ) is the characteristic exponent of a 
subordinate Brownian motion. Since ip(X) « X a ' 2 A A, we have if; G WUSC(a, 1, C) n 
WLSC(«, l,c) for some < c 1 s= C. 

7. Let < ax < a 2 < 2 and u(r) = r" 1 / 2 " 1 Vr^ 2 ' 1 . Let 77(A) = Cu{X) = A^ l/2 7(ai/2, A) + 
A- a2 / 2 r(a 2 /2,_A) and <p(\) = l/rj(\). Note that <p(X) « X ai ^ A A a2//2 . Therefore <^(|x| 2 ) G 
WUSC(a 2 , 0, C) and WLSC(«i, 0, c) for some < c ^ 1 < C. It is shown in [261 Example 
10.18(i)] that if is a special Bernstein function but not a complete Bernstein function. 
Moreover, the Levy measure of if is not known and so previous methods of estimating 
transition densities of the resulting subordinate Brownian motion do not yet apply [IS] . 
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4.2 Estimates 

The following estimate is a version of [33j Theorem 7 (ii) (b)] with explicit constants. 

Lemma 13. Let f ^ be nonincreasing, (3 > and Cf G WUSC(— 0, 9, C) . There is b 
b(/3,C) G (0,1) such that 

f(r) ^ ^r^Cfir- 1 ), < r < b/9. 
Proof. Let < b < 1. If u > 9, then by Lemma [5] and the upper scaling (with A = /u), 

pbu^ 1 poo 

uCf(u) = u e- us f(s)ds + u e- us f(s)ds 

JO Jbu" 1 
pbu~ 



< " 



7 (1 



fbu 1 poo 

L — / e- us £f(s- 1 )s~ 1 ds + fibu- 1 ) / e- us uds 

i 1) JO Jbu- 1 



pbu 



tTTI) j ° ^ US ^ C f( u ) e ~ USs ~ 1(ls + f{bu- l )e- b 

= C^\uCf(u) + f(bu- l )e- b . 
7(1,1) 

If 2C^(f3,b) ^ 7(1, 1) = 1 — e -1 , then fibu^ 1 ) ^ e b uCf(u)/2. We change variables: r = ira -1 . 
Since Cf is decreasing, 

/(r) ^ ^eV" 1 /:/^^ 1 ) ^ ^-^/(r- 1 ), r < b/9. 

□ 

Lemma 14. C = C(d) exists such that if ipE WUSC(a,9,C) and a = [(2 — a)C] 2 , then 

F(\X t \ ^ r) > a (1 - e -^* (1/r) ) , < r < 

Proof. By Lemma HI Proposition [2] and ( JT4"j) . for A ^ 1 and 9 we have, 



rfJXv) 1 _ p -7T 2 tV(VA^) 1 _ -7r 2 #(v^)CA Q / 2 

f^^Z <; ^A' 1 1 — ^ < C^ 1 — ^ vr 2 ^ 2 ^ 2 - 1 . 

£/t(u) 1 l-e-^v^) ^ 1 i_ e -*^(V5) ^ 1 



Thus, Cf t GWUSC(a/2 - 1, 9 2 , n 2 CfC). By Lemma US and Lemma H 

¥(\X t \ >r) = f t (r 2 ) > b -e b r- 2 Cf t {r~ 2 ) > (l - e -***(Vr)) f r ? < h/ e\ 

Here 6 G (0, 1) is such that 27r 2 C 2 C7(l — a/2, b) ^ 1 — e -1 , see the proof of Lemma [131 Since 
7(1 - a/2, b) < b l - a ' 2 /{l - a/2), we may take 6 = ( ^^=J and a = 6/(2Ci) < 1. □ 

Since lim t _ >0 + ^(1^1 ^ r )A = = L(r) for r > 0, we obtain the following result. 

Corollary 15. //^ satisfies WUSC(a,9,C) and a is from Lemma 14, then 

L{r) ^ aip*{l/r), <r < y/E/6. 
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We recall that a reverse inequality is valid more universally, cf. ffTTj) . 

In view of f|T8|) and (120]) . the lower weak scaling implies power- type asymptotic growth of 
ip and so we can use Fourier inversion to estimate p t . Here is a preparation. 

Lemma 16. Let a > and \1/ be an increasing function on [0, oo), such that ^(O) = and 
^ EWLSC(a,6,c). There is C 2 = C 2 {d,a) such that if < K < oo, t > and tV(8) < K, 
then 

[ e- m ^d£ ^ C 2 (1 V (cK)-^- 1 ) (*- 1 (K/t)) d . 
Proof. Since \l/ increases and scales, it is unbounded. Hence, for t > 0, 

„ oo „ oo 

JR d n=1 J cK(n-l)^{\t\)<cKn « n=1 

Also, if #(0) < i*T/t, then ^/~ l (K/t) > 6. By lower scaling, for n ^ 1, 

^(n 1 /"^ 1 ^/^) ^ = cnK/t, 



which yields 
We obtain, 



n 



V a %-\K/t) > y-\cKn/t). 



J 



d 

n=l 



For p ^ and u > we have 



u{n+l) 



n=l n=l 
/■oo 

< e 2u / x p e- ua: rfx = e 2u M" 1 ^r(p+ 1). 



Since S(u,p) is decreasing in w, we also have S(u,p) ^ 5'(l,p) for w ^ 1, thus 

5(«,p)^e 2 r(p + l)(lV«" w ), «>o. 
The proof is complete: J2n=i n^e"^™" 1 ) ^ e 2 r(rf/a + 1)(1 V (cK)-*/"' 1 ) . □ 

For a continuous nondecreasing function : [0, oo) — >■ [0, oo), such that 0(0) = 0, we let 
0(oo) = lim^oo cj)(s) and we define the generalized inverse 0~ : [0, oo) — > [0, oo), 

4>~{u) = inf{s ^ : 0(s) ^ w}, ^ u < oo, 

with the convention that inf = oo. The function is nondecreasing and caglad where finite. 
Notice that <p(<p~(u)) = u for u G [0, 0(oo)] and _ (0(s)) ^ s for s G [0,oo). Also, if <p : 
[0, oo) —7- [0,oo), <p(0) = 0, c > and c0 ^ <p, then (j)~(u) ^ <p~(cu), u ^ 0. Below we shall 
often consider the (unbounded) characteristic exponent ip of a unimodal Levy process with 
infinite Levy measure and its (comparable) maximal function x/j*, and denote 

r = (rr- 

This short notation is motivated by the following equality, 

inf{s ^ : ^{s) ^ u} = inf{s ^ : ip*{s) ^u}, ^ u < oo. 

Note that ip* (u)) = u, ip~ (ip* (s)) ^ s. The reader may find it instructive to prove the 
following result for t > and x G lR d \ {0}. 
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Lemma 17. tip*{l/\x\) > 1 if and only if tip*{l/\x\)\x\- d ^ [if>~ (1 / 't)) d . 

In what follows it may be helpful to view ttp*(l/\x\) ^ 1 and tip*(l/\x\) < 1 as defining 
"large time" and "small time" for given x G M. d \ {0}, respectively. 

We note that ip~ (tp* (s)+) ^ s for s G [0, oo), where ?/>~(u+) denotes the right hand side 
limit of at u ^ 0. Furthermore, scaling of ?/> translates into scaling of ip~ as follows. 

Lemma 18. Iftp G WLSC(a, 0, c) ; tfienV~ G WLSC (l/2, 0, (c/t 4 ) 1 /^ W(l/a, 0, (vr 3 /c) 2 ^) . 

Proo/. We let *(s) = s > 0. Lemma HU and Remark Q] yield * G WUSC(2,0,1). By 

Lemma [T]and Remark [2], \I/ G WLSC(a, 0, c/tt 2 ). Remark [3] implies 

tf -1 G WLSC(l/2, 0, 1) n WUSC(l/a, 0, (n 2 /c) l/ ^). 

By Lemma [1] and the above scaling, 

* _1 (s/2) ^ ^~(s) < ^- 1 (stt 2 /2) < (tt 4 /c) 1/ ^~ 1 (s/2). 

Hence, by Remark |2] G WLSC(l/2, 0, (c/tt 4 ) 1 ^) n WUSC(l/a, 0, (tt 3 /c) 2 ^). □ 

We shall use Fourier inversion and ([T]) to estimate Pt(0): if, say, lim^^^ ■?/;(£)/ ln(|£|) = oo, 
then e~*^® is integrable for t > 0, M. d 3 x h> ^(a;) may be assumed continuous, and we have, 

Pt (x) = (2n)- d [ e-^e-^df. 

In particular, 

pt(0) ^ (2vr)- d / e-^* (l5l) ^ > ^y d — [i/;-(l/t)] d , t > 0. (22) 

Lower scaling yields a reverse inequality. 

Proposition 19. If a > and ip G WXS'C(a,^, c), then C = C(d,a) exists such that 

Pt(x) ^ Cc-^ - 1 [ip-(l/t)] d if t > and t^*{6) < 1/tt 2 . (23) 

Proof. We let 

*(s) = ^(s" 1 ), s>0. 
Note that /ii and \l/ are strictly monotone. According to Lemma [U, 

2^(s)/tt 2 ^ tp(s) <: 4>*(s) 2*(s), 

hence \l> G WLSC(a, c/7r 2 ). Furthermore, i/;~(2u/ir 2 ) ^ vf- 1 ^) ^ ij)-{2u) for w ^ 0. Let 
t > 0. If ty>*(0) < lA 2 , then (2t/7T 2 )^(^) < 1/tt 2 . We apply Lemma HI to X = 1/tt 2 and 
(time) 2t/vr 2 : 

p t (x) < (27r)- d / e-^Ciei)/^^ < (27r)- d C 2 (c/vr 2 )^ 1 [^(^i)- 1 )]" . 
But ^ _1 ((2t) _1 ) ^ ip~(l/t), and the proof is complete. □ 
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Under the assumptions of Proposition [T9l by (l22|) and (1231) we obtain 
e- 1 < ft(0) : [r [l$* Wd < e 2 r(d/a + 1) (vr 2 /c) d/a+1 , * > 0, *0*(g) < 1/tt 2 . (24) 

(We may interchange pt(0) and [^~(l/t)] d in formulas.) 

Also, if a > and ip G WLSC(a, 0, c), then for t > 0, x E R d , 

p t (x) ^ C(d,a)c- d / a - 1 [ip~(l/t)] d . (25) 

We note in passing that the same argument covers the Gaussian case = |£| 2 and more 
general exponents otherwise excluded from our general considerations. We also note that 
analogues of ff25|) are often obtained by using Nash inequalities [3j [2] . 

For the a-stable process addressed in (j4]), by using directly Lemma [T6l with K — 1 we get 

2e 2 r(^) 
Pt[X) 2 d K d / 2 aT{d/2) 

and the constant is not too far off from the exact estimate, obtained by direct integration; 



Pt{x) ^ pt(0) = (2tt) 



-d 



I* /*00 

/ e -^ a di = (2n)- d / w^^-^dr 

JR d JO 



oo 



1 2r (-) 

= (27r)- d uj d I s d l a - l -t- dta e- s ds = . ^f- r r d/a . (26) 
K J J a 2 d n d / 2 aT(d/2) v ; 

Note that in this case ip~(u) = u x l a . The above change of variables tr a = s nicely counterpoints 
the proof of Lemma [T6l and the role of the quantity tip*(l/\x\) appearing, e.g., in Lemma [T71 

Corollary 20. If ij) eWLSC(a,9,c) with 9 > 0, and < T < oo ; then p t (x) < C[ip-(l/t)) d 
for allO <t <T and x e R d , with C = C(d,T,a,9,c,i/;). 

Proof. Note that ip* satisfies WLSC(a, 9, c/tc 2 ). Since ip*(0) = and ip* is continuous and 
unbounded, there is 9 1 > such that ip*(9 1 ) = l/(n 2 T). By Lemma[l2l ip* satisfies WLSC(a, #A 
9i,C\), with a constant c\. An application of Proposition [T9l completes the proof. In fact, if 
ip e WLSC(a,0,c), then C = C{d,T,a,9,c,ip{9)). Indeed, if 9 1 < 9, then l/(vr 2 T) =ij)*{B x ) sC 
^(0){c/tt 2 )- 1 {9/9 1 )-^, hence 0f > c^-/[vr 4 T^(^)], leading to C = C(T, a,0,c,V>(0)). □ 

Thus, (I2"5j) holds for all i > 0, even if 9 > 0, but the constant deteriorates for large t. 

The following main result of our paper gives common bounds for unimodal convolution 
semigroups with scaling. Notably, our second main result, Theorem [26] below, shows in addition 
that scaling is equivalent to common bounds. 

Theorem 21. If ip eWLSC(a,9,c), then there is C* = C*(d,a) such that 

Pt(x) ^ C*min jfe)-^- 1 [i)-(l/t)] d , ^^p^ l if * > and tip* (9) < 1/vr 2 . 
IfipE WLSC(a, 9,c)nWUSC(a, 9, C), then c* = c*(d,a,c,a, C), r = r (d,a,c,a, C) exist with 
Pt(x) > c* min J fyr (1/t)] d , |VJ } */ *>0, tip*(9/r ) < 1 and |z| < r o /0. 
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Proof. Let t > 0. For x = the term in the statement should be ignored-the 

bounds are to be understood as f [22l and fl23|) . Accordingly, below we let i 6 1^ {0}. The 
upper bound now follows from Corollary [7J and Proposition [TUJ 

To prove the lower bound, we consider k ^ 2, to be determined later. We have 

\B K \ x \\B lxl \ u d (K d -l) 

Let |x| < y/a/6, with a = a(d,a,C) from Lemma [T41 We now suppose that tip*(l/\x\) ^ 1. By 
concavity, — e _s ^s for ^ s ^ 1. By Lemma [T41 and Corollary El 

¥(\X t \^\x\)-F(\X t \^K\x\) > >(l/|x|)-^-(2 ( i+ l)^*(l/|«x|) 

2 e — 1 

^ 4e(2d + l)V*(l/M) 



2 T w 1 " V a(e-l) 

Recall that ^* GWLSC(a, 0, c/vr 2 ). If |x| < 1/(k0), then 

4e(2d + 1) V*(VM) < 47r 2 e(2d +1) _ 1 
a(e-l) ^ ca(e - 1) K ~ 2' 

where 

« = (8n 2 e(2d + l)/(ca(e - 1))) V - (27) 

Recall that c, a 6 (0, 1] (see the proof of Lemma [HI), so k ^ 2, as required. Also, k~ x < ^[a. 
Therefore this choice of k yields 

Pt(x) > , j - jtr{V\x\)\x\- d (28) 
(jJd{K a — I) 4 



4:U d (K d - 1) [' W ^ ' \x\ d 

We are in a position to verify that the lower bound in the statement of the theorem holds 
with r = minjfir 1 , y/a} = kT 1 . We thus assume that t > 0, tip*{6/r ) < 1, and because of 
the preceding discussion we only need to resolve the case < \x\ < r /6, tip*(l/\x\) ^ 1. By 
continuity, there is x* e R d such that |x| < \x*\ < r /9 and tip*(l/\x*\) = 1. By ( 128]) . 

Pt(x) > Pt (x*) > -^A? > 7^7 



ad . /r / _ /l/ .x 1 d tip*{l/\x\) \ 



This ends the proof and we may take c* = adn /{Audj- □ 



Remark 5. For the record we note that the constants in the lower bound depend on d, a, c, a, C 
via ( 1271) and Lemma fT41 e.g. c* = c(d, a, a) c d ^- 2 ^ d+g - ) ^ 2 -^ anc [ rQ _ c (a,a,d) (c 2 cj 

In view of Lemma [T71 the two factors in the minima in the statement of Theorem [211 should 
be interpreted as the approximations of Pt{x) in large time (on-diagonal regime) and small time 
(off-diagonal regime), correspondingly. 

We emphasize that the upper bound in Theorem [21] only requires the lower scaling. For 
instance the upper bound holds for the 2-regularly varying characteristic exponent = 
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|£| 2 /[log(l + I'Cl 2 )]^ with G (0,1), which is in agreement with the outcome of the Davies' 
method in this case |22j . 

Tracking constants for the isotropic a-stable Levy process addressed in (J3J), we obtain 



Pt{x) > c a ( t~ d l a A ; 



\x 



d-} 



where 

2(d+a) „ j , _, , v (4-a)(ri+oQ 

(2-«)«(l-e- 1 )^r(f,i)^^ 



q , 3d , (8-a)(d+a) d 2d (12-2a)(d+t») d . 2(d+a) , 

2 ~ <x(2-a) 7] -2+~ H c«(2-a) (2d + 1) a + <*(2-<*) F 1 f|j a(2-a) 

We now list a number of general consequences of Theorem EH We first complement ffl6|) by a 
similar lower bound resulting from Theorem 1211 

Corollary 22. 1/^ G WLSC(a, 9, c)D WUSC(a, 9, C) and\x\ < r /9, thenu(x) ^ c*4)*{l/\x\)\x\- d . 
Corollary 23. If ip £ WLSC(a,0,c)nWUSC(a,0,C), then (E} ZioWs /or aW t > and x G R '. 
Proof. We use ( TToT) and Corollary [22] to obtain 



» " ' y , iGr\{0}. (30) 

We then appeal to (122]) . (1231) and Theorem EH □ 

By scaling, in particular by Lemma [TBI we obtain the following important doubling property, 
cf. [T7] in this connection. 



Corollary 24. Ifip satisfies (global) WLSC(a,0,c) and WUSC(a,0,C), then 
Pt(2x) ~ Pt(x) and P2t{ x ) ~ Pt{ x )i t > 0, x G R d . 

Thus, if 6* = in Theorem I2T1 then the global asymptotics of pt(x) is fully and conveniently 
reflected by tp. If 9 > 0, then our bounds are only guaranteed to hold in bounded time and 
space (bounded time for the upper bound) . For large times we merely offer the following simple 
exercise of monotonicity. 

Corollary 25. IfipE WLSC(a, 6_, c) , < \x\ < (tt- 4 c) V - /9 and tip* ({x^ 1 ) ^ 1, then 

pt(0) ^ C(d,a)c- d fe- 1 til)*(\x\- 1 )/\x\ i . 

Proof. Define (threshold time) t = l/ip*^" 1 ). By Proposition El ip* GWLSC(a, 9, c/tt 2 ), 
thus 

Since t —> Pt(0) is decreasing, by Proposition liTHwe have for t ^ to, 

Pt(0) < Pto (0) < Cc-^- 1 ^-(r(\x\- 1 ))] d < Cc- d ^- l -\- d ^ Cc~ d l 

which completes the proof. □ 
The next theorem proves that our definitions quite capture the subject of the study. 



tip*(\x 


n 




X ' 


i 
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Theorem 26. Let X t be an isotropic unimodal Levy process in IR d with transition density p, 
Levy-Khintchine exponent ip and Levy measure density v. The following are equivalent: 

(i) WLSC and WUSC [global WLSC and WUSC] hold for ip. 

(ii) For some r G (0, oo) [r = oo] and a constant c, 

pAx) ^ c y )\ 1 ; , < Id < r Q , < tif>*(\x\ _1 ) < I. 
\x\ d 

(Hi) For some r G (0, oo) [r = oo] and a constant c, 

, N ^(kr 1 ) 

v(x) > c , ' , < x < r . 



Proof. Theorem I2T1 and Lemma [T7l yield the implication (z) =>■ (ii). The implication (ii) =>- (iii) 
follows because lim t _ s>0 + x )/t — v { x ) vaguely on IR d \ {0}. To prove that (iii) implies (i), we 
assume that (iii) holds. By the Levy-Khintchine formula, if)(x) = cr|x| 2 + f Rd (l— cos(£, x))z/(<ix), 
where a ^ 0. Actually, we must have cr = 0, because 

f |2 /" .nV'd^l^ 1 ) i /" C i 

oo > / |x| v(x)dx ^ / c|x| — j— — ax ^ a ——^dx. 

JBx JB 1ArQ \ x \ JB 1ArQ \ x \ 

By [261 proof of Theorem 6.2] and (JSJ), the following defines a complete Bernstein function, 



<p{\) 



poo \ /*oo 

/ -^s- 1 v{s- 1/2 )s- d/2 ds= / (1 - e- Xu ) n(u)du, A^0, 



where n(u) = C[v(s l l 2 )s d l 2 }{u). In fact, by changing variables, and (jH]) for A > we have 

poo \ 2 /-oo 

</?(A) = 2 / t/(u)u <t - 1 dM « / [1 A (Am 2 )] viuju^du = u^h (A _1/2 ) . 

Jo Am + 1 Jq 

By Corollary [HI , there exists c\ = c\(d) such that 

Ciip(X) ^ if) (y/\j ^ q V(A), A ^ 0. (31) 
Since <p'(\) = J Q °° ue~ Xu n(u)du and /i is decreasing, Lemma [5] with n = and m = 1 yields 

7(2,1) u 

Using the upper incomplete gamma function and monotonicity of u, we obtain 

e- s \ x \\{s' l / 2 )s- d / 2 ds w (| x |2) 
" (x) ^ * jj~ a e --W V*/ 2 da ^ |x|^ 2 r(l-rf/2,l)' X ^°- (33) 

We leave it at that for a moment, to make another observation. As ip is a complete Bernstein 
function, we have that 

^(A) = AMA) 
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is a special Bernstein function (see [26| Definition 10.1 and Proposition 7.1]). Since X t is pure- 
jump, limi^ooVKO/KI 2 = 0. Thus, lin^^ <fi(\) = oo. Also, <p(0) = 0, and by [261 (10.9) 
and Theorem 10.3], the potential measure of the subordinator with the Laplace exponent <pi is 
absolutely continuous with the density function, 

/i(u)du. 

In particular Cf = l/(pi. Let x £ R d be such that < \x\ < r Q . By (jHTJ, (Hi), fl33|) and fl32|) . 



ClC ^£!l < c t^l < „ (a0 < MM!) < 5ffl3 ( 34) 

1 " \x\ d " w " |x| d - 2 r(i - d/2, i) ^ |x| rf + 2 r(i - d/2, i) 7 (2, i) ■ { v 

Therefore C2<p(\) ^ Ay/(A) for A > l/^/ro, where c 2 = ccij(2, 1)T(1 — d/2, 1). This implies 
that the function A~ C2 <^(A) is nondecreasing on (l/y/r~o, oo). By Remark [TT1 and ( 13TI) . </? G 
WLSC(c 2 , l/v^d, 1)- Hence, V GWLSC(2c 2 , r^ 1 , Cl ). 

We now prove the upper scaling. By concavity of ip, u<p'(u) ^ <p(u). For < s < ^Jr\ i1 by 



/(s) > c 3 y <^(m x )u x tiw ^ c 3 y <p'(« l )u 2 du = c 3 (<p(s x ) - y?(r 1/2 )) . 

Note that ip is strictly increasing. Therefore, for < s < ^ro/2, we get 

f{s) > c^is- 1 ) = a/lstp^s- 1 )]. (35) 
Since / is decreasing and Cf(u) = l/ipi(u), by Lemma [5] we obtain 



7 (2, l)s 2 \<pj " ' 7 (2,l)s 2 ^( S - 



Thus, by (|35l) . C5(pi(X) ^ Ay/^A), where A > 2/y / ro and C5 = 047(2, 1). It follows as above 
that (p x e WLSC(c 5 ,2r - 1/2 , 1). Since <pi is concave, \<pi(\) ^ ^i(A), hence C5 < 1. Thus, 
up E WUSC(1 -c 5 ,2r - 1/2 ,l). This implies ip e WUSC(2(1 - c 5 ), ^ \ cf). □ 

The reader may consult, e.g., example 3 in Section I47T1 for typical asymptotics of pt(x) and 
v(x). Our next observation results from Corollary [22], ( 1301) . ( 13TI) and Theorem [261 

Corollary 27. If the characteristic exponent of a unimodal (isotropic) Levy process X satisfies 
global WLSC and WUSC (with exponents < a ^ a < 2), then there is a complete subordinate 
Brownian motion with comparable characteristic exponent, Levy measure and transition density. 

To indicate an application of Corollary [271 we remark that the boundary Harnack principle 
proved in [Tj5] for complete subordinate Brownian motions under global scaling conditions 
should now extend to general subordinate Brownian motions with global scaling. 

We see from Corollary [27J that the asymptotics of the characteristic exponent, Levy measure 
and transition density observed for complete subordinate Brownian motions are representative 
among all unimodal Levy processes with lower and upper global scaling. 

We close our discussion with a related result, which negotiates the asymptotics of the Levy 
density (at zero) and the Levy-Khintchine exponent (at infinity) under approximate unimodality 
and weak local scaling conditions. We hope the result will help extend the common bounds. 

Recall our conventions: 0<a,57<2, 0<c^l^C<oo. 
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Proposition 28. Let X be a pure-jump symmetric Levy process with Levy measure u(dx) = 
u(x)dx and characteristic exponent ip. Suppose that 9 G [0, oo), constant c G (0, 1] and nonde- 
creasing function f : (0, oo) — y (0, oo) exist such that 

/(VM) , x i/(VN) 
c^LU ^ u(x) ^ c" 1 ^^, < |x| < 1/9. 

Iff G WLSC(a,9,c)n WUSC(a,9,C), then /(|£|) andip(£) are comparable for |f | > 9. In fact, 
there is a complete subordinate Brownian motion whose characteristic exponent is comparable 
to ip(x) for |£| > 9, and whose Levy measure is comparable to v on Bi/g \ {0}. 

Proof. Let Y be the pure-jump unimodal Levy process with Levy density 

v ¥ (x) = f(l/\x\)\x\- d l Bl/e (x) } xeR d \ {0}. 

Let ip Y be the characteristic exponent of Y. Using C = C(d) of ( IT6|) and Proposition [21 we get 

/deix^c^e), teR d \{o}. (36) 

On the other hand Corollary [3] yields C = C'(d) such that 



^ Y (O^C'[ [(|£lM) 2 Al]^Mcfe, ^G 
J B 1/B \ z \ 

By scaling of /, for |£| > 9 we have 

^ y (0 <c7(iei) (Wi 2 - 3 / n ^5 +e " 1 ^"° f dz 



^l^ + AWi)- (37) 



2 — a ca / 
By the symmetry u(z) = v(—z), 

^(0=/ (1-005(^,^)^)^-^(0+/ (l-cos^M^tfc, (38) 



(Si/s) c 

in particular ^ w ?/; y on M d if = 0. If 9 > 0, then the last integral in ( 13 8 p is bounded by 
2v ((i?i/6i) c ) < oo; by Proposition [2] we have, 

Vtt) > ^P Y (0)/n 2 > 0, \x\ > 9, 

and so ip(£) ~ t/> y (£) ~ /(|£|) f° r \C\ > ®, where in the second comparison we used fl36l) and 
d37D. It follows that i) Y G WLSC(o, 9, c t ) n WUSC(o, 0, (7 X ). By (ETJ) in the proof of Theorem 
1261 there is a complete subordinate Brownian motion Z whose characteristic exponent ip z is 
comparable with ip Y on IR d . Its Levy density v z is comparable with z/ y on B ro /g \ {0} by ( fTBT) 
and Corollary [22l where ro = ro(e£, a, c 1; a, Cx). The comparability of v z [x) and v Y [x) also 
takes place on \ B ro / e because the functions are bounded from above and below on the set, 
as follows from (T2T]) and monotonicity of / > 0. Thus, v and v z are comparable on B\iq. □ 

To clarify, the semigroup of the process X in Proposition [28] is not necessarily unimodal, 
hence its estimates by / call for other methods, e.g. those based on u, mentioned in Section [TJ 
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